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We theoretically analyze the eflFect of electron-electron interactions on Aharonov-Bohm (AB) cur- 
rent oscillations in ring-shaped systems with metallic quantum dots pierced by external magnetic 
field. We demonstrate that electron-electron interactions suppress the amplitude of AB oscillations 
Iab at all temperatures down to T = and formulate quantitative predictions which can be ver- 
ified in future experiments. We argue that the main physical reason for such interaction-induced 
suppression of Iab is electron dephasing while Coulomb blockade effects remain insignificant in the 
case of metallic quantum dots considered here. We also emphasize a direct relation between our 
results and the so-called P(i5)-theory describing tunneling of interacting electrons. 



PACS numbers: 72.10.-d 

I. INTRODUCTION 

Aharonov-Bohm (AB) oscillations of conductance as a 
function of the magnetic flux 4> piercing the system rep- 
resent one of the fundamental properties of meso- and 
nanoscale conductors which is directly related to quan- 
tum coherence of electronsi. Coherent electrons propa- 
gating along different paths in multiply connected con- 
ductors, such as, e.g., metallic rings, can interfere. Such 
interference effect results in a specific quantum contri- 
bution to the system conductance 5G. Threading the 
ring by an external magnetic flux $ one can control the 
relative phase of the wave functions of interfering elec- 
trons, thus changing the magnitude of 5G as a function 
of The dependence (5G($) turns out to be periodic 
with the fundamental period equal to the flux quantum 
$0 = hc/e. 

It is important to emphasize that the phase of the elec- 
tron wave function is sensitive to its particular path. In 
diffusive conductors electrons can propagate along very 
many different paths, hence picking up different phases. 
Averaging over these (random) phases or, equivalently, 
over disorder configurations yields the amplitude of AB 
oscillations 5G{^) with the period 'I'o to vanish in diffu- 
sive conductors-^. There exists, however, a special class 
of electron trajectories which interference is not sensi- 
tive to disorder averaging. These are all pairs of time- 
reversed paths which are also responsible for the phe- 
nomenon of weak localization^. In multiply connected 
disordered conductors interference between these trajec- 
tories gives rise to non- vanishing AB oscillations with the 
principle period $o/2- Such oscillations will be analyzed 
below in this paper. 

It is well known that various kinds of interactions. 



such as electron-electron and electron-phonon interac- 
tions, electron scattering on magnetic impurities etc. can 
lead to decoherence of electrons thus reducing their abil- 
ity to interfere. Accordingly, AB oscillations should be 
sensitive to all these processes and can be used as a tool 
to probe the fundamental effect of interactions on quan- 
tum coherence of electrons in nanoscale conductors. Re- 
cently it was demonstrated^i^i^ that the effect of quan- 
tum decoherence by electron-electron interactions can be 
conveniently studied employing the model of a system 
of coupled quantum dots (or scatterers). This model 
might embrace essentially all types of disordered conduc- 
tors and allows for a straightforward non-perturbative 
treatment of electron-electron interactions. It also allows 
to establish a direct and transparent relation^'^ between 
the problem of quantum decoherence by electron-electron 
interactions and the so-called P{E) theory^ii, see also^ 
for an earlier discussion of this important point. In this 
paper we employ a similar model in order to study the 
effect of electron-electron interactions on AB oscillations 
in disordered nanorings. 

The structure of our paper is as follows. In Sec. 2 
we define our model and outline our general real time 
path integral formalism employed in this work. Sec. 3 
is devoted to a detailed derivation of the effective ac- 
tion for our problem in terms of fluctuating Hubbard- 
Stratonovich fields mediating electron-electron interac- 
tions. With the aid of this effective action we then evalu- 
ate Aharonov-Bohm conductance of the ring in the pres- 
ence of electron-electron interactions. This task is ac- 
complished in Sec. 4. A brief discussion of our results is 
presented in Sec. 5. Some technical details of disorder 
averaging are relegated to Appendices. 
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FIG. 1: The ring-shaped quantum dot structure under con- 
sideration. 



II. THE MODEL AND BASIC FORMALISM 

Below we will analyze the system depicted in Fig. 1. 
The structure consists of two chaotic quantum dots (L 
and R) characterized by mean level spacing Sl and Sr. 
Here we will restrict our attention to the case of metallic 
quantum dots with Sl,r being the lowest energy param- 
eters in the problem. These dots are interconnected via 
two tunnel junctions Ji and J2 with conductances Gti 
and Gt2 forming a ring-shaped configuration as shown in 
Fig. 1. The left and right dots are also connected to 
the leads (LL and RL) respectively via the barriers Jl 
and Jfl with conductances Gl and Gr. We also define 
the corresponding dimensionless conductances of all four 
barriers as (741,2 = Gti,2Rq and gL,R = Gti,2Rq, where 
Rq = 27r/e^ is the quantum resistance unit. These dimen- 
sionless conductances are related to the barrier channel 
transmissions Tk via the standard formula g = 2^j,Tfc, 
where the sum is taken over all conducting channels in 
the corresponding barrier and an extra factor 2 accounts 
for the electron spin. 

For the sake of convenience in what follows we will 
assume that dimensionless conductances gL,R are much 
larger than unity, while the conductances gti and gt2 are 
small as compared to those of the outer barriers, i.e. 



gL,gR > 1,. gti, 3*2- 



(1) 



The whole structure is pierced by the magnetic flux $ 
through the hole between two central barriers in such way 
that electrons passing from left to right through different 
junctions acquire different geometric phases. Applying a 
voltage across the system one induces the current which 
shows AB oscillations with changing the external flux $. 

The system depicted in Fig. 1 is described by the ef- 
fective Hamiltonian: 



H 



i,j=L,R 

+ Hj+Tl+Tr + T, 

j=L,R 



(2) 



where Cy is the capacitance matrix, V£(/f) is the electric 
potential operator on the left (right) quantum dot. 



H 



LL 



d'r^^rLirmLL-eVLL)'^,, .LL(r), 



Hrl=Y. J d''r^l,RLir){HRL~eVRL)^c.Mr) 

are the Hamiltonians of the left and right leads, Vll,rl 
are the electric potentials of the leads fixed by the exter- 
nal voltage source, 

defines the Hamiltonians of the left {j = L) and right 
(j — R) quantum dots and 



2m 



^l + U,{r) 



is the one-particle Hamiltonian of electron in j-th quan- 
tum dot with disorder potential Uj(r). Electron transfer 
between the left and the right quantum dots will be de- 
scribed by the Hamiltonian 



E 



J1+./2 



Here the integration runs over the total area of both tun- 
nel barriers Ji and J2. The Hamiltonian Tl{r) describing 
electron transfer between the left dot and the left lead 
(the right dot and the right lead) is defined analogously 
and are omitted here. 

Before we proceed with our analysis the following two 
remarks are in order. Firstly, we point out that within 
our approach the effect of electron-electron interactions 
is accounted for by the voltage operators Vl,r in the 
effective Hamiltonian In order to avoid misunder- 
standings we would like to emphasize that this approach 
is fully equivalent to one employing the usual Coulomb 
interaction term in the initial Hamiltonian of the system. 
The operators Vl.r. corresponding to fluctuating poten- 
tials of the left and right dots emerge as a result of the 
exact Hubbard-Stratonovich decoupling of the Coulomb 
term containing the product of four electron operators. 
This is a standard procedure (described in details, e.g., 
in Ref. and elsewhere) which is bypassed here for the 
sake of brevity. 

Secondly, we note that in Ref. we have studied weak 
localization effects in a system of coupled quantum dots 
within the framework of the scattering matrix formalism 
combined with the non-linear tr-model. However, in order 
to incorporate interaction effects into our consideration 
- similarly to Refs. 4,5 - it will be convenient for us to 
describe inter-dot electron transfer within the tunneling 
Hamiltonian approach, as specified above. For clarity let 
us briefly recapitulate the relation between these two ap- 
proaches. For this purpose we define the matrix elements 
tim = mTlm) between the l—th wave function in the left 
dot and m— th wave function in the right dot. Electron 
transfer between these dots can then be described by a 
set of eigenvalues of this matrix tk where, as above, the 
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index k labels the conducting channels. These eigenval- 
ues are related to the barrier channel transmissions 



Q 

as* 



Tfc 



(3) 



This equation allows to keep track of the relation between 
two approaches at every stage of our calculation. 

We now proceed employing the path integral Keldysh 
technique. The time evolution of the density matrix of 
our system is described by the standard equation 



(4) 



where H is given by Eq. Let us express the opera- 
tors e^'^* and e''^* via path integrals over the fluctuating 
electric potentials V^'^ defined respectively on the for- 
ward and backward parts of the Keldysh contour: 



~iHt 



iHt 



DVf T exp <^ -i / dt'H [Vf{t')] 



= DVfTc^plt dt'H[Vfit')]\. (5) 



Here T exp (T exp) stands for the time ordered (anti- 
ordered) exponent and the Hamiltonians H [V^^ (i')] , 

[^^(t')] ^re obtained from the original Hamiltonian 
^ if one replaces the operators Vj(t) respectively by 
the fluctuating voltages Vf{t') and Vf{t'). 

Let us define the effective action of our system 



iS[V^,V^] = In tr 



Texp<^-i / dt'H[vf{t')] 



X /5oTexp<^i / dt'H[vf{t')] 



(6) 



Since the operators H [Vf'{t')] , H [Vj^{t')] are quadratic 
in the electron creation and annihilation operators, it is 
possible to integrate out the fermionic variables and to 
rewrite the action in the form 



iS = iSc + iSext + 2Trln [G^^] 



(7) 



Here Sc is the standard term describing charging ef- 
fects, Sext accounts for an external circuit and G""*" is 
the inverse Green-Keldysh function of electrons, moving 
in fluctuating voltages field. It has the following matrix 
structure: 



G 



(Gil Tl 

ft 





f 



\ 



T; 



R 






Tr 
Gr\ I 



(8) 



Here each quantum dot as well as two leads is repre- 
sented by the 2x2 matrix in the Keldysh space: 



idt -H, + eV/ 






-idt + m 



(9) 



Tunneling blocks has the following structure in Keldysh 
space: 



Tl,r — 







Jl,r 



J tL,R{r')5{r' -r)dr' 
(10) 







T 



' J t{r')5{r' ~r)dr' 
J1+J2 

/ t{r')6{r' - r)dr' 

J1+J2 



III. EFFECTIVE ACTION 



(11) 



In what follows it will be convenient for us to remove 
the fluctuating voltage variables and the vector potential 
from the bare Green functions. This is achieved by per- 
forming a unitary transformation under the trace in Eq. 
([7]). As a result we find 



(12) 



-e'v^ J t{r')S{r' - r)dr' 







T, 



^i^B J t{r')S{r' - r)dr' 
Ji 



Here we introduced the fluctuating phase differences 



(13) 



'PF,B{t) 



dT{Vl^{r)^V['''{r)) (14) 



defined on the forward and backward parts of the 
Keldysh contour as well as the geometric phases 



c?x^j4^t(a;). 



(15) 



where the integration contour starts in the left dot, 
crosses the first (ifii^^) or the second (^'g^'') junction and 
ends in the right dot. The difference between these two 
geometric phases equals to 



(1) 



= 2^$/$o, 



(16) 



where $ is the magnetic flux threading our system. 

Let us now expand the exact action iS ^ in powers 
of T. Keeping the terms up to the fourth order in the 
tunneling amplitude, we obtain 



iS « iSc + iSext + iSL + iSr - 2tr 



GlTGrT^ 



tr GLTGRf^GLTGRf^ 



..(17) 





FIG. 2: Diagrammatic representation of different contributions originating from expansion of the effective action in powers of 
the central barrier transmissions: second order (AES) terms (a) and different fourth order terms (b,c). 



The terms iSl.r define the contributions of the isolated 
dots (which are of no interest for us here), the second or- 
der terms oc yield the well known Ambegaokar-Eckern- 
Schon (AES) action^ iS"^^^ , and the fourth order terms 
oc i'* account for the weak localization correction to the 
system conductance-*^. 

Let us first analyze the AES action. Performing aver- 
aging of this action over disorder in each dot separately 



as well as averaging of tunneling amplitudes with the 
correlation function 



we arrive at the following result 



iS^^' ^ - j dt,dt2 J rfx^-fj^ Yl G)^(xti;xt2)(-l)^e^'^^(*^)G]l(xi2;xti)(-l)'e-^-(*^), (19) 

where the convention (—1)^ = —1, (—1)^ = 1 is implied. This AES contribution to the action is described by the 
standard diagram depicted in Fig. 2a. We observe that after disorder averaging the AES action becomes totally 
independent of the magnetic flux. Hence, this part of the action does not account for the AB effect investigated here. 

In order to evaluate the contribution sensitive to the magnetic fiux $ it is necessary to analyze the last term in Eq. 
p?)l . Averaging over realizations of transmission amplitudes yields two types of terms illustrated by the diagrams 
in Fig. [2)d,c. It is straightforward to check that only the contribution generated by the diagram (c) depends on the 
external magnetic fiux, while the diagram (b) does not depend on $. On top of that, the terms originating from the 
diagram (b) turn out to be parametrically small for metallic quantum dots considered here. This observation will be 
justified in Appendix A. 

It follows from the above arguments that only the diagram in Fig. 2c is responsible for the AB effect in our system. 
Its contribution to the action reads 

f J. J. J. f [ ^u,9t(^)9t(y) ^ 



iS^ = - g2z(>p(")->p(")) j dtidt2dhdU j dx j dy 



X J2 G'2'(xti;y<2)(-l)^e'^^-(*^)Gi,'=(yt2;xt3)(-l)'e-*'^'=(*^) X 

i,jMJ=F,B 

xG|'(xt3;yt4)(-l)'e''^'(*^^G^,(yt4;xii)(-l)V-*'^'(*i). (20) 

Since Gl,r are the equilibrium Green-Keldysh functions of the dots they can be expressed via retarded (G^) and 
advanced {G^) Green functions in the standard manner: 

GL,fl(xiii;x2i2) = J rft(Gf^fl,(xifi;x2t)Fi(t-<2)- A(ti -i)Gf_^(xii;x2^2)), (21) 



where 



Here /(<) = / f{E)dE/ 271 is the Fourier transform of the Fermi function f{E) ~ (exp(i?/T) -1-1) ^ and h{t) = 

5{t)-m. 
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What remains is to combine Eqs. ([21]) and (|20|) and to average the latter over disorder. This procedure amounts 
to evaluating the averages of the products of retarded and advanced Green functions in each dot separately. Such 
averaging can be conveniently accomplished either by means of the diagram technique or with the aid of the non- linear 
CT- model. The corresponding calculation is presented in Appendix A. It yields (i — L, R): 

(Gf (xiti;X2t2)Gf (X3t3;x4t4))d = (Gf (xiii;X2i2))<i(Gf (X3t3;x4i4))d, (23) 



(Gf (Xiti; X2t2)Gf (Xgtg; X4t4))d (Gf (xi^i ; X2i2))d(Gf (xgts; X4t4))d, 



(24) 



(Gf(xiii;X2t2)Gf (X3t3;x4i4))d = (Gf(Xiti;X2i2))d(Gf (X3t3;x4t4))d + 

+27rA'"in7(|xi — X4|)ti7(|x2 — xsj) x 

Xi + X4 X2 + X3 ~ 



6{ti-t2 + h-U) + 



xCi ti - t2] 



2 ' 2 
+27r7V,n7(|xi - X3|)tu(|x2 - X4I) x 

Xi + X3 X2 + X4 



(25) 



where I?L^i?(t; x, y) and CL,fl(t;x, y) the diffusons and the Cooperons in the left and right dots and vj{r) = 
e"*"/^' sin fci?r//ci?r. Substituting these averages into the action (|20p it is straightforward to observe that only the 
terms containing the product of two Cooperons yield the contribution which depends on the magnetic flux $. This 
part of the action takes the form 



ra.n— 1.2 



2»(^r'-^r') / dr,dr2 / dtMt.dU I dK I dy^^fe^ x 



Jr, Jrr 



xC,(ri; y, x)CH(r2; x, y)e^(.+ (*.)-^+(*3)+.+ ^in "LM x 

, / N ..- y~ ('2) „, , (*2) 

h{ti ^t2~ Ti)e'^- + /(ti - i2 - Ti)e-'^- 



h{t2 ~h- T2)e~'^^f{t3 -u + n)- 



-f{t2 - <3 - T2)e''^^h{h ~u + n) 



e'^^fiU ~h+ T2) + e~'^^h{U -ti+ T2) 



(26) 



where we defined the "classical" and the "quantum" com- 
ponents of the fluctuating phase: 



^+(0 - 



V3f(0 + •^B{t) 



LP (t) = ipF{t) - ipsit). 

(27) 

The above expression for the action S^'" (pS)) fully ac- 
counts for coherent oscillations of the system conduc- 
tance in the lowest non-vanishing order in tunneling. It 
is important to emphasize that no additional approxima- 
tions were employed during its derivation and, in par- 
ticular, the fluctuating phases are exactly accounted for. 
We will make use of this fact in the next section while 
considering the effect of electron-electron interactions on 
AB oscillations in the system under consideration. 



IV. CURRENT OSCILLATIONS 



Let us now evaluate the current / through our system. 
For this purpose we will employ a general formula 



2,„± 



6S[ip+,ip ] rS\v+.v-] 



(28) 



Substituting the total effective action into this formula we 
arrive at the result for the current which can be split into 
two terms / = /q + 51 ^ where /g is the flux-independent 
contribution and 51 is the quantum correction to the cur- 
rent sensitive to the magnetic flux (f>. This correction is 
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determined by the action iS^^, i.e. 



61 



(29) 



Below we will only be interested in finding the quantum 
correction (^5)1 . 

In order to evaluate the path integral over the phases 
ip^ in (|29p we note that the contributions Sc and 5cxt in 
Eq. ([T7|) are quadratic in the fluctuating phases provided 
our external circuit consists of linear elements. Other 
contributions to the action are, strictly speaking, non- 
Gaussian. However, in the interesting for us here metal- 
lic limit HI) phase fluctuations can be considered small 
down to exponentially low energiesi^iil in which case it 
sufflces to expand both contributions up to the second or- 
der (/?^. Moreover, this Gaussian approximation becomes 
exac^ in the limit of fully open left and right barriers 
with gL,R ^ 1. Thus, in the metallic limit ([1]) the inte- 
gral l|29p remains Gaussian at all relevant energies and 
can easily be performed. 

This task can be accomplished with the aid of the fol- 
lowing correlation functions 



(ifi+it)) ^eVt, (^-(t))-0, 



(30) 



((^+(t)-^+(0))^+(0))=-i^(t), (31) 



{v+{t)^-{0) + ^-{tW+m = 2iK{\t\), (32) 



(^+(i)^-(0) - v-{t)v^m = 2zX(i), (33) 



(^-(t)^-(O)) =0, 



(34) 



where the last relation follows directly from the causality 
principle^^. Here and below we define V = Vrl — Vll to 
be the transport voltage across our system. 

Substituting the AB action ([^5]) into Eq. (|29p one ar- 
rives at the expression containing six different phase av- 
erages listed in Appendix B. All these averages in Eqs. 



(IBip - (IB6p are expressed in terms of two real correla- 
tion functions F{t) = {{^{t) - <^(0))2)/2 and K{t) = 
i{[(f{0),<f{t)])/2 defined above in Eqs. §1^ and ([5^ . 
Note that these correlation functions are well familiar 
from the so-called P(i?)-theory®'^ describing electron 
tunneling in the presence of an external environment 
which can also mimic electron-electron interactions in 
metallic conductors. They are expressed in terms of an 
effective impedance Z{uj) "seen" by the central barriers 
Ji and J2 



— coth— K[Z(w)] (35) 



m^e^J'^nZi.)f-^. (36) 

Further evaluation of these correlation functions for our 
system is straightforward and yields 



Fit) 



In 



sinh(7rrt) 



27r 

K{t) ~ — sign(i), 



(37) 



(38) 



where we defined g = 47r/e^Z(0) and 7 ~ 0.577 is 
the Euler constant. Neglecting the contribution of ex- 
ternal leads (which can be trivially restored if needed) 
and making use of the inequality ([1]) we obtain g ~ 

25i5i?/(3L +5i?)- 

We observe that while F{t) grows with time at any 
temperature including T = 0, the function K(t) always 
remains small in the limit g 3> 1 considered here. As 
we demonstrate in Appendix B, the correlation function 
F{t) should be fully kept in the exponent in Eqs. (|B1[) - 
()B6P while the correlator K{t) can be safely ignored in 
the leading order in 1/g. Then combining all terms we 
observe that the Fermi function f{E) - though present 
in the effective action (|26|) - drops out from the final ex- 
pression for the quantum correction to the current which 
takes the form: 



m.n—L.R 

xCi(ri;y,x)CH(r2;x,y)e-2^(^^)-2^(^^)+^(^-^^)+^(-^+^^). (39) 



We observe that the amplitude of AB oscillations is af- 
fected by the electron-electron interaction only via the 
correlation functions for the "classical" component of the 
Hubbard-Stratonovich phase ip'^ . Both the correlators 



containing the "quantum" phase ip~ and the Fermi func- 
tion f{E) enter only in the next order in 1/g which de- 
fines weak Coulomb correction to 61 ignored here. For 
more details on this point we refer the reader to Refs. 
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m 

The result ([55]) can also be rewritten as 

«($) = -lABm ~ IWLI - IWL2, (40) 

where the first - flux dependent - term in the right-hand 
side explicitly accounts for AB oscillations and reads 

Iab{'^) = Iab cos(47r$/$o), (41) 

while the last two terms IwLi,2 represent the remaining 
part of the quantum correction to the current which does 
not depend on $. 

Already at this stage we would like to clarify the rela- 
tion between our present results for AB oscillations and 
those for WL correction to conductance^. In order to 
derive Eq. ((39)) we have evaluated the contributions of 
all processes illustrated by the diagrams in Fig. 2b, c 
and identified terms sensitive to the magnetic field which 
were not considered in Ref. Q. In this way we have ob- 
tained the AB current lAsi'^) in Eqs- HOI), dH]) which 



represents our new result to be analyzed below. The 
two remaining terms in Eg. (HO)) are the WL corrections 
already evaluated in Ref. 0. Towards the end of this 
section we will explicitly specify the relation between all 
three contributions to the quantum correction ()40p . 

Let us evaluate the amplitude of AB oscillations Iab 
for the system with two identical quantum dots with vol- 
ume V, dwell time td and dimensionless conductances 
9l — gR = g — 47r/(5T£), where 5 ~ 1/VN is the dot mean 
level spacing and N is the electron density of states. In 
this case the Cooperons take the form 



Ci(i;x,y) =C;i(t;x,y) = -^e-"^- . (42) 



Defining dimensionless conductances of central barriers 
as gti,2 = 2 gt{x)dx we obtain 



Iab — 



e'^gtigt2S'^V 
47r3 



dTidT2e 



-2F{ti}-2F(t2)+F{ti-T2)+F{ti+T2) 



(43) 



In the absence of electron-electron interactions {F{t) 
0) this formula yields: 



AO) 
^ab 



4e^gtigt2^ 



(44) 



In order to account for the effect of interactions we 
need to specify the effective impedance Zioj). Its real 
part takes the form 



47r / T 
^Z[u) = — [- 



1 



1 



7rS{uj) 
td + trc 



(45) 



where 1/t = I/td + I/trc, trc — T^/gEc is the i?C-time 
and Ec is an effective charging energy of our system. Eq. 
(|43)) demonstrates that electron-electron interactions al- 
ways tend to suppress the amplitude of AB oscillations 
Iab below its non- interacting value (|44p . Combining 
Eqs. ([37)) and (gS)) with ([43]) at high enough temper- 
atures we obtain 



Iab 



-^ {2^TTncf'<' 



•D 



T\") I 1 ( gTRC \ ^/'^ 

^AB K 2td y T ) ' 



Tu^ < T. 



.-1 

' RC 



(46) 



while in the low temperature limit we find 

^AB \ u / 



•d 



(47) 



The latter result demonstrates that interaction-induced 
suppression of AB oscillations in metallic dots with 
'Trc ^ 'Td persists down to T = 0. 

The ratio Iab / Iab evaluated numerically as 

a function of temperature at different values of g. The 
corresponding results are presented in Fig. [3l We observe 
that - in accordance with the above analytic expressions 

- the ratio Iab / Iab grows with decreasing T as a power 
law and finally saturates to a constant value smaller than 
unity at r < 1/td. The suppression of AB oscillations 

- both at higher temperatures and at T ^ clearly 
depends on the interaction strength which is controlled 
by the parameter 1/g in our model. Fig. [4] demonstrates 
the dependence of Iab / ^ab 9 limit of zero 
temperature and for td/tjic — 10. While at moderate 
values oi g 10 20 interaction-induced suppression 
of Iab remains pronounced down to T = 0, at weaker 
interactions {g > 100) this effect becomes less significant 
and is merely important at higher temperatures, cf. Fig. 

H 

In order to complete our analysis let us briefly address 
additional quantum corrections to the current IwLi,2 in 
Eq. (j^O)) . Although these terms do not depend on <& 
and, hence, are irrelevant for AB oscillations, they allow 
to establish a direct and transparent relation between the 
Aharonov-Bohm effect studied here and the phenomenon 
of weak localization in systems of metallic quantum dots 
with electron-electron interactions'''^. With the aid of 
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DQ 
< 




FIG. 3; (Color online) The ratio Iab/i'"2b versus temperature 
at different values of dimensionless conductance g. 



m 
< 




FIG. 4: (Color online) The ratio Iab/Iab 

as a function of 
dimensionless conductance g a.t T = and td/trc ~ 10. 



Eq. ((39|) one easily finds 



IwLi 
Iab 



9ti 
25*2' 



IwL2 

Iab 



9t2 



(48) 



Combining this equation with the above resuhs for Iab 
we immediately identify the terms IwLi and Iwl2 as 
weak localization corrections to the current^'^ originat- 
ing from the two central barriers in our structure. In 
addition, in the absence of the magnetic field $ = the 
total quantum correction to the current (5/(0) (|40)) ex- 
actly coincides with the weak localization correction to 
the current for two connected in series metallic quantum 
dots^!^ provided the two central barriers in Fig. 1 are 



viewed as a composite tunnel barrier with total dimen- 
sionless conductance gn + gt2- 



V. CONCLUDING REMARKS 

The established relation between our present results 
and those obtained in Refs. 00 helps to clarify the main 
physical reason for the effect of interaction-induced sup- 
pression of AB oscillations in our structure. In full anal- 
ogy with the weak localization correction^ both at non- 
zero temperatures and T = this suppression is due 
to electron dephasing by electron-electron interactions. 
This decoherence effect reduces the electron ability to 
interfere and, hence, decreases the amplitude Iab below 
its non-interacting value /^"^ . At the same time Coulomb 
blockade effect - although yields an additional suppres- 
sion of Iab ~ remains weak in metallic quantum dots and 
can be neglected as compared to the dominating effect 
of electron dephasing. It is also important to empha- 
size that in the course of our analysis we employed only 
one significant approximation: We performed a regular 
expansion of the current in powers of the tunneling con- 
ductances up to second order terms (forth order terms in 
the tunneling matrix elements). At the same time the ef- 
fect of electron-electron interactions on AB oscillations in 
our system was treated non-perturbatively to all orders 
and essentially exactly. 

Note that one could be tempted to interpret the sup- 
pression of Iab aXT — Q just as a result of a simple renor- 
malization effect by electron-electron interactions which 
is not related to dephasing. It is important to stress that 
- unlike, e.g., in the case of the interaction correction for 
single quantum dot o^^i^^ - here such interpretation would 
not be appropriate. The fundamental reason is that the 
interaction of an electron with an effective environment 
(produced by other electrons) effectively breaks down the 
time-reversal symmetry and, hence, causes both dissi- 
pation and dephasing for interacting electrons down to 
T = Oi^. In this respect it is also important to point out 
a deep relation between interaction-induced electron de- 
coherence and the P(S)-theory^ which we already em- 
phasized elsewhere^'^'* and which is also evident from our 
present results. Similarly tc'''^ one can also introduce the 
electron dephasing time in our problem and demonstrate 
that at r — > it saturates to a finite value in agreement 
with available experimental observations ^^'^^i^^ . We be- 
lieve that the quantum dot rings considered here can be 
directly used for further experimental investigations of 
quantum coherence of interacting electrons in nanoscale 
conductors at low temperatures. We also note that our 
model can possibly be applied to analyze the behavior of 
recently fabricated self-assembled quantum rings^^ where 
the AB oscillations have been observed by means of mag- 
netization experiments. 
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FIG. 5: Diagrammatic representation for vertices r£j(xi, X2; w), rc(xi,X2;a') and for averages Xci(xi,X2;£) and Xt;(xi, X2; e). 
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APPENDIX A: AVERAGING OVER DISORDER 

Let us consider the following disorder averages of the product for retarded and advanced Green functions for one 
of the quantum dots: 

Xd(xi,X2;£) = (G^(xi,X2;w)G'^(x2,xi;cj - e))d - (G^(xi, X2; cj))d(G'^(x2, xi; cj - e))<j, (Al) 
Xc(xi,x2;e) = (G^(xi,x2;a;)G^(xi,x2;a;-e))d - (G^(xi, X2; w))d(G'^(xi, X2; w - £))d, (A2) 
Xfc(xi,X2;e) = (G^(xi,xi;tj)G^(x2,x2;t^-e))d- (G^(xi,xi;t^)>d(G'^(x2,x2;w-£))d, (A3) 



where 



G^(-^)(xitr;x2t2) = /^e— (*-*^)G^(-^)(xr,X2;^). 



(A4) 



In order to evaluate the above averages we will employ the standard diagram technique for noninteracting electrons 
in disordered systems'^'^. The essential elements here are the so-called diffuson and Cooperon ladders depicted in Fig. 
[5] where we also define vertices rd(xi,X2; w) and rc(xi,X2; w). In the presence of time-reversal symmetry and in the 
limit of low momenta and frequencies these vertices obey a diffusion-like equation: 



-ia;-i:»V^jr<j(c)(xi,X2;tj) = 2ttNt'^ ^^^^ ~ 



(A5) 



Here D = vpl/i and = l/vp are respectively the diffusion coefficient and the electron elastic mean free time. With 
the aid of the above vertices one can define the diffuson and the Cooperon respectively as 



/did 
^e-^"*rd(xi,x2;w), 



(A6) 



C(i;xi,X2) = 27r7VT| / ^e-^"'re(xi, X2; 
J 27r 



CO). 



(A7) 



In the absence of the magnetic field they obey the following diffusion equations 



{dt - DVlmt; xi, X2) = (5(xi - X2)(5(t), 



(A8) 



{dt ~ DS/lJC{t; xi, X2) = (5(xi - X2)<5(0 (A9) 
with appropriate boundary conditions. 
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FIG. 6: Diagrams which define the average Xfe(xi, X2; e). 



Evaluating the diagrams for X(;(xi,X2;e) depicted in Fig. [5] after some algebra we arrive at the following result 

Xi + X2 Xi + X2 



where 



Xd(xi,X2;e) = (27r7Vre)"rd(xi,X2;e) + (2^iVTe)^n7^(|xi -X2|)r, 



tn(|xi -X2I) = J dx(G^(xi,x;w))d(G^(x,X2;w))d. 



(AlO) 



(All) 



In the case of 3d systems we find zu{r) = e sinkpr /kpr. 
The expression for Xc(xi,X2;£) is derived analogously. We find 

X,(xi, X2; e) - (27rA^Te)2r,(xi, X2; £) + {2^NT,fw^{\^^ - X2|)rd 



/Xl + X2 Xl + X2 

V 2 ' 2 



(A12) 



Combining Eqs. jJT]), (El (EH, JMI, jM} with (|XT0l) - (|XT2)) we arrive at Eq. (^5]) . 

Note that the average Xfc(xi,X2;e) (|A3[) is omitted in Eq. ([25]) since this average turns out to be parametrically 
small as compared to both Xj^{'x.^, X2; e) and A"c(xi, X2; e). In order to demonstrate this fact it is necessary to evaluate 
the diagrams for Xfe(xi,X2;£) depicted in Fig. [6l Proceeding as above we get 



Xfc(xi,X2;e) {2TlNTe) W (|xi -X2 



/xi + X2 Xl + X2 



Xl + X2 Xl + X2 



2 ' 2 

{2TiNTlf [r2(xi,x2;e) + r2(xi,X2;e)] . (A13) 



The first term in this equation clearly vanishes for |xi — X2I > /. Here we assume that the size of both the dots and 
the contacts is large as compared to the electron mean free path I. Provided the typical contact size is of the same 
order as that of the dots the latter condition implies I <C vftq- If, however, the contact size is much smaller than 
that of the dots this condition becomes ppi ^ V^ch , where Nch is the effective number of conducting channels in 
the contact. In this case both the diffuson and the Cooperon do not depend on coordinates and are defined by Eq. 
(li^ . Then one finds 

Substituting this result into Eq. (|XT3l) we get Xk{t) w V^^9{t)te-*^/''° . Comparing this expression with that for 
Xd{c){t) « 27rNV-'^e{t)e-^^''° at times t < td we obtain 



Xk/Xa(c) < td/{NV) ^ td5 ^ l/g « 1. 

This estimate demonstrates that the average Xk (|A3[) can be safely disregarded in Eq. 
consideration. 



(A14) 

for the problem under 



APPENDIX B: AVERAGING OVER FLUCTUATING PHASES 



Substituting the action (|26p into Eq. (P^j) one expresses the current as a combination of different phase 

averages evaluated with the total action As we already argued above, in the metallic limit (H]) all these 
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averages are essentially Gaussian and, hence, can be easily performed. For the sake of completeness, below we present 
the corresponding results: 



^~F(ti-t2)-F{ti-t4.)~F{t2-t3)-F(t3-t4)+F{ti-t3) + F{t2~t4) ^ 
^^-iK{ti-t2)-iK{\ti-t3\)-iK{ti-U)+iK{t2-t3)+iKi\t2-U\)-iK(t3-t4) 

^-F{tl-t2)-F{tl-U)-Fit2-t3)-F(t3-t4)+F{tl-t3)+F{t2-ti) ^ 

^^~iK(ti-t2)-iK{\ti-t3\)-iK{\ti-ti\)+iK{t2-t3}-iK{t2-ti)+tKi\t3-ti\) ^g2) 

^g»(V+(t2)-V+fe) + V+(t4)-V'+(tl) + ^^ + ^^-^^ + ^^)) = 
^-F{ti-t2)-F{ti~U)-F{t2-t3)-F(t3-t4)+F{ti-t3)+F{t2-t4) ^ 
y^^-iK{ti-t2)+iK{ti-t3)-lK(ti-U)-iK{\t2-t3\) + iK{\t2-t4\)-lK(\t3-t4\) (^23-) 

^gZ(V+(t2)-¥'+fe)+'P+(t4)-V+(tl) + ±^S^-^^T^ + ^^-^^)) = 
^-F{ti^t2)-F{ti-t4)-F{t2-t3)-F{t3-U) + F{ti-t3) + F{t2-ti) ^ 

y.^-iK{ti-t2)-iK(\ti-t3\)+iK{\ti-t4\)+iK{\t2-t3\)-iK{\t2-t4\)+iK{\t3-U\) ^g^^ 

^g«(v>+(t2)->P+(t3) + ¥'+(t4)-¥'+(tl) + ^^S^ + ^'^T^-^^T^ + ^^^)) = 
g-F(ii-t2)-F(ti-t4)--F(t2-t3)--F(t3-i4) + F(ti-t3) + F(i2-i4) X 
y^g-iK{\tl-t2\)+iK{\tl-t3\)-iK(\tl-t4\)-iK(\t2-t3\)+iK(\t2-t4\)-iK{\t3-t4\) ^gg^ 

^g,:(^+(t2)-v'+(t3)+v+(t4)-v+(ti)+^^ + ^^-^^-^^)) = 

g-F(ii-t2)-F(ti-t4)-F(t2-i3)-F(t3-t4)+F(ti-t3) + F(t2-t4) y. 
y.^-iK{ti-t2)+iK{ti-t3)+lK(\ti-t4\)+iK{\t2-t3\)-iK{\t2-t4\) + iK{\t3-t4\) |-gg^ 



Note that for arbitrary metallic conductors gL.R. 3> 1 
all these equations are accurate down to exponentially 
small energies (set by the so-called renormalized charg- 
ing energyiSiiii which is of little importance for us here) , 
and in the particular limit of fully open left and right 
barriers Eqs. (jBip - (jB6P become exact^^. Thus, combin- 
ing (|Bi])-(Il6| with Eqs. ([26]), dMl), (137]) and (138]) we 
exactly account for the effect of electron-electron inter- 
actions on the amplitude of AB oscillations in the system 
under consideration. 

It is useful to observe that in order to quantitatively 
describe this effect in the metallic limit g ^ I one can to- 



tally neglect aU the functions K{t) in all Eqs. (|BT1) - (|B6|1 . 
This is because these functions remain much smaller than 
one at all times (cf. Eq. ([55]) ) and, hence, can only cause 
a weak (~ Coulomb correction to Iab which fur- 

ther slightly decreases the amplitude of AB oscillations. 
The origin of this Coulomb correction is exactly the same 
as that identified and discussed in the weak localization 
problem^. Thus, no additional discussion of this point 
is necessary here. 

Substituting unity instead of all the exponents in (jBip - 
(|B6|) containing iiT-functions and keeping all F-functions 
in the exponent, one easily arrives at Eq. (j39p . 
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